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Abstract 

We measure the frequencies of small-amplitude shape oscillations of a magnetically-levitated 
water droplet. The droplet levitates in a magnetogravitational potential trap. The restoring forces 
of the trap, acting on the droplet's surface in addition to the surface tension, increase the frequency 
of the oscillations. We derive the eigenfrequencies of the normal mode vibrations of a spherical 
droplet in the trap and compare them with our experimental measurements. We also consider the 
effect of the shape of the potential trap on the eigenfrequencies. 



'Electronic address: 



richard . hill@nottingham.ac . uk 



I. INTRODUCTION 



If the surface of a spherical liquid drop is briefly deformed, by a puff of air, for example, 
it vibrates, ringing at several different frequencies. The eigenfrequencies of these shape 
oscillations were determined by Lord Rayleigh, 



where T is the surface tension, p is the density and a is the radius of the spherical drop at 
rest [TJ |2] . By measuring a, we can determine the surface tension of the liquid. Beaugnon 
et al. used diamagnetic levitation to measure a of the lowest order (7 = 2) mode of a 
diamagnetically levitated liquid droplet [3] and, recently, we used diamagnetic levitation to 
investigate dynamics of a spinning water droplet |4J. A diamagnetically- levitated droplet 
is confined within a magnetogravitational potential trap [6] . The trap acts as an addi- 
tional cohesive force on the drop, perturbing its eigenfrequency spectrum. Beaugnon et al. 
observed the shift to higher frequency of the lowest-order I = 2 mode, and determined an 
expression for the increase in terms of an enhanced effective surface tension [3]. However, 
their result cannot be generalized to the higher order modes. Here, we use diamagnetic 
levitation to measure the normal mode frequencies of a levitating drop, for / > 2. We derive 
an expression for the eigenfrequencies of a liquid droplet confined by the magnetogravita- 
tional potential trap and compare it with our measured frequencies. Our analysis points the 
way to achieving accurate measurements of surface tension using this non-contact technique. 
We also consider the effect on the eigenfrequencies of the shape of the potential trap and 
consider analogies with the vibrations of a model 'star'. 

II. EXPERIMENTAL DETAILS 

We use a vertical-bore superconducting solenoid magnet with a room-temperature, 50 mm 
diameter bore to levitate droplets of water with radii ~ 1 cm. The droplets levitate approx- 
imately 80 mm above the geometric center of the solenoid, where the diamagnetic force, 
proportional to BVB, is equal in magnitude to, and opposite in direction to, the gravita- 
tional force on (i.e. the weight of) the droplet (5j E]. The magnetic field is B « 12 T and 
the vertical field gradient is dB/dz f» 120 Tm" 1 at the levitation point. 




(1) 



2 




-20 -15 -10 -5 5 10 15 20 
x (mm) 



FIG. 1: Bottom: Magnetogravitational potential U(x,z) of a water droplet levitated by the mag- 
netic field B(x, z) of a vertical-bore superconducting solenoid magnet (contours at AU / g = 0.05 mm 
intervals); x and z are radial and vertical cylindrical coordinates, respectively, with origin at the 
geometric center of the solenoid. The field at the center of the solenoid is Bq = B(0, 0) = 16.5 T. 
The field profile B(x, z), used to construct the plot, was computed by numerical integration of the 
Biot Savart integral. The stable levitation point at a local minimum in the potential is labeled Lq. 
An unstable levitation point at a saddle-point in the potential is labeled L\. Top: schematic of the 
experimental set-up. The magnetic field lines are shown in gray. Rays are shown as bold lines. 



Fig. [T] shows a spatial map of the magnetogravitational potential 
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of a unit mass of water in the potential trap, where x and z are radial and vertical cylindrical 
coordinates with origin at the geometric center of the solenoid coil, x = ~ 9 x 10~ 6 and 
p = 1 x 10 3 kgm -3 are the SI volume magnetic susceptibility and density of water respectively, 



and g = 9.8 ms~ 2 . The center of the levitating droplet coincides with the position of the 
stable levitation point L at the center of the trap, as shown. (There is a second point of 
unstable levitation L\, at the saddle point in U, as shown in Fig. [TJ However, this point is 
not useful for these experiments, and so we do not discuss it further here.) The net force on 
a unit mass of water, T = —VU, is zero at the levitation point. The equilibrium shape of 
a liquid with no surface tension follows the contours of U, as demonstrated recently, using 
liquid H 2 close to the critical point [7] . For water droplets with a ~ 1 cm, however, the 
surface tension dominates the magnetic and gravitational forces on the drop, so that its 
equilibrium shape is nearly spherical. 

The shape of the magnetogravitational potential trap can be altered by adjusting the 
current in the magnet solenoid coils. A convenient measure of the current is the magnetic 
field at the geometric center of the solenoid, Bq = B(0,0). The field B(x,z) is everywhere 
proportional to B Q . To levitate water requires B Q ps 16 — 17 T using our magnet. Fig.[T]shows 
the U(x, z) for B = 16.5 T. The field profile B(x, z) of the magnet, used to construct this 
plot, was computed by numerical integration of the Biot Savart integral, using a thin-shell 
approximation for the current density in the solenoid. 

We expand the potential U in a multipole expansion about Lq\ 

U(r,0) = Y,Ci(r)P j (coBe), (3) 
i>o 

where r and 9 are spherical coordinates with origin at Lo; 9 is the polar angle (i.e. rsin# = 
x). Only the derivative of U normal to the droplet's surface (i.e. the radial component of 
the force) influences the eigenfrequencies of the normal modes, 

IVM) = = J2c>(r)P 3 (cos9), (4) 

j>o 

where c'- = dcj/dr. By adjusting B , we reduce the quadrupole component until it is small 
compared to the spherically-symmetric component (c' 2 <C c' ). The octopole harmonic (c' 3 ), 
which cannot be reduced this way, remains comparable to c' . All other harmonics are small 
compared to c' . We shall return to discuss these points further in Sec. |V} 

Fig. [l] top, shows a schematic diagram of the water droplet levitating in the vertical 
magnet bore. The droplet at rest is close to spherical; we measure the ratio of the equatorial 
(horizontal) diameter to the polar (vertical) diameter to be 1.00 ± 0.02. An optical fiber 
directs light from a HeNe laser at the droplet. The drop focuses the light onto the aperture 
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of a second fiber, which transmits the light to a photodiode outside the magnet. The EMF of 
the photodiode is measured by a storage oscilloscope. A 1 mm-diameter nozzle, directed at 
the center of the underside of the droplet (see Fig. 1), is connected to a rubber bulb outside 
the magnet by a tube. When the bulb is struck on a hard surface, the resulting pulse of 
air from the nozzle excites several shape oscillation modes simultaneously, with amplitude 
< 0.05a. Since the focal length of the drop depends on its shape, the intensity of the laser 
light falling on the photodiode oscillates as the drop vibrates. The temperature of the water 
was brought to 16°C in a water bath before it was injected into the potential trap, to match 
the ambient temperature in the magnet bore. At this temperature, the density of the liquid 
is p = 999 kgm -3 , its surface tension is T = 73.3 mNm -1 and its kinematic viscosity is 
v = 1.11 x 10 -6 m 2 s _1 [H]. The liquid was injected into the trap using a glass pipette. The 
volume of liquid injected was determined to better than 1% uncertainty from the difference 
in the weight of the pipette before and after injection. The liquid was drawn out of the bore 
using a paper towel, by capillary action, after the experiment. By measuring the difference 
in weight between the wet and dry paper, we obtained a second measurement of the droplet 
volume. Using this simple and accurate technique, we were able to determine that the 
mass loss through evaporation of the drop, during the measurement period (approximately 
30-60 minutes), was always less than 2%. Experiments were performed both in air and 
in dry nitrogen gas (by filling the bore with N 2 from a pressurized gas cylinder). The air 
experiments were performed at B = 16.2 T and the nitrogen experiments at 16.5 T. The 
reason for the different B may be explained by the oxygen content of the air, which, being 
paramagnetic, buoys up the droplet (9J [TO] by an additional force r air = — VC/ a i r . Here, 
£4h- = Xair-B 2 /(2yU p) and x a i r = +3.7 x 10~ 7 is the volume magnetic susceptibility of air 
(S.I. units). The c' 2 component of U + t/ a i r is minimized at Bq = 16.2 T. In nitrogen, the 
buoyancy force is negligible compared to the diamagnetic force, and c' 2 is minimized at 16.5 T 
as shown in Fig. [T] 



III. RESULTS 



Fig. [2J bottom, shows the oscillations in the photodiode EMF £(t) developed by the light 
refracted through an a = 6 mm droplet; they decay exponentially with a time constant r ~ 
0.1 — 10 s, dependent on a and I, due to the viscosity of the water, i.e. £{t) = J2i £i exp(ic^t — 
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FIG. 2: Bottom: photodiode voltage £ measuring the shape oscillations of an a = 5.8 mm water 
droplet. Center: power spectrum of the oscillations of an a = 5.8 mm water droplet, as a function 
of oscillation frequency / = u/2tt. Top: measured frequencies / of the oscillation modes of water 
droplets with radii between a = 4.5 mm and 14 mm, up to I = 16. The broken lines show the 
Rayleigh frequency spectrum for a water drop, according to Eqn. [TJ The measured frequencies 
(crosses) are slightly higher, by ~ 1 Hz, due to the effect of the magnetogravitational potential 
trap. 

t/ri). The power spectrum of these oscillations is shown in Fig. |2j center. Several peaks 
are evident in this spectrum, corresponding to the I — 2 to 7 Rayleigh modes of the drop. 
Peak 1 is due to a small oscillation of the droplet's center of mass of about L Q (amplitude 
~ 1 mm). The frequency of this peak is independent of a. The upper panel of Fig. [2] shows 
the measured peak positions for different drops with radii between a = 4.5 mm and 14 mm; 
the broken lines show the Rayleigh frequency spectrum for a water drop, according to Eqn. [TJ 
Note that the measured frequencies are slightly higher, by ~ 1 Hz, due to the effect of the 
magnetogravitational potential trap. We shall calculate this frequency shift in the following 



sections of the paper. Although small, it is important to be able to account for this shift, if 
the technique is to be used to obtain accurate measurements of surface tension, for example. 

An estimate of the effect of viscosity on the frequency of oscillation can be obtained 
by calculating the ratio of the magnitude of viscous stresses to surface tension forces, 
fvp/iTD- 1 ) = Oh[2/(/ - 1)(Z + 2)] 1 / 2 / 7 r, where D — 2a, f — 2iru and Oh is the Ohne- 
sorge number [11]. This ratio is much smaller than unity for the droplet sizes used in these 
experiments. This indicates that we can neglect the influence of viscosity on the oscillation 
frequencies as a small effect. For example, we expect the viscosity of the water to marginally 
lower the frequency u of the / = 2 mode of an a = 5 mm droplet by (5 ± 1) x 1CT 4 %, but 
this is a small reduction compared to the frequency increase resulting from the trapping 
potential. (We obtained these estimates by solving numerically the Chandrasekhar equa- 
tion for the eigenfrequencies of a viscous drop [I2| 113]). The Oh number indicates that 
viscosity has a significant effect on the frequency of modes / < 20 only for water drops 
smaller than a ~ 10 /mi. The shape of the peaks in the power spectrum agrees well with 
the Lorentzian shape expected for an exponentially decaying oscillation, with half-width 
(HWHM) Aw = 1/r ^ va- 2 (l - \){2l + 1) given by the Chandrasekhar equation [121 113]. 
(Since the power spectrum is the square of the magnitude of the Fourier transform of the 
oscillations, we compare the shape of the peak with the square of the Lorentzian function). 



IV. SPHERICAL POTENTIAL APPROXIMATION 

We now consider the effect of a spherically-symmetric magnetogravitational potential well 
on the eigenfrequencies of the drop. We shall discuss the effects of additional harmonics in 
Sec. |Vj Our derivation follows Lamb's derivation of the Rayleigh frequencies [2] closely, with 
the addition of the force T r on the droplet 's surface due to the gradient of the magnetograv- 
itational potential at the surface. We write the shape of the Ith mode of a drop oscillating 
with frequency u, for oscillations with small amplitude e, as 

r = R(6,t) = a + ( = a + ePi(cos6) smut (5) 

where Pi is a Legendre polynomial of degree / > 1 (/ = 1 corresponds to an oscillation of 
the droplet's center of mass in the potential trap). 
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The pressure equation at the surface of the droplet (to first order) is [2] 

d(j> 
di 



U{R) + ^- + F(t), (6) 



where F(t) is an arbitrary function of t only and is the velocity potential [2] 



<p(r,6) = — - j—£ uePi(cos 9) cosut. (7) 
The pressure difference across the surface resulting from the surface tension is [2] 



a a 

The magnetogravitational potential U at the surface of the drop is (to first order) 

U(R(9, t)) = U(a) - T r (a, 6)(R - a) (9a) 
= U(a) - eT r (a,9)Pi(cos9) sia(ut), (9b) 

where T r (a, 9) = —cf (a) for a spherically symmetric well (see Eqn. [5]). Inserting Eqns. |9bJ [8] 
and [7] into Eqn. [6j we obtain 

to 2 = 4 + a 2 , (10) 

where 

2 c (a)Z 

= ^ ' 

is the oscillation frequency (squared) of a hypothetical drop with T = 0, held together by 

the magnetogravitational trap alone. The fact that the square of the measured frequency 

uj 2 is a simple sum of the square of the Rayleigh frequency a\ and a 2 is due to the fact 

that I remains a good eigennumber for oscillations in a spherically- symmetric potential well. 

In Sec. [V] we will consider the effect of a non-spherical magnetogravitational well. In a 

non-spherical well, I is not a good eigennumber in general, but for small deviations from 



spherical, we can use perturbation theory to obtain corrections to Eqn. 10 

We can obtain an experimental measurement of c' by examining the difference between 
the measured frequencies u of any two of the modes I, n > 2: 

where hi = /(/ — 1)(/ + 2). Dividing by (nhi — lh n )/a, we obtain an experimental measurement 
of c' (r) from the oscillations of a drop that has radius a = r at rest: 

"4fr)=r^^. d3) 
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FIG. 3: Left panel: Gradient of the magneto-gravitational potential trap c' obtained from the 
measured frequencies of two modes, I = 2 and I > 2, of a drop with radius a = r at rest. Error 
bars show the standard error. Filled circles: data obtained in nitrogen atmosphere. Open squares: 
data obtained in air atmosphere. Dotted line: c computed from the solenoid geometry. Right 
panel: dependence of the experimentally-obtained Cq on the mode number with I > 2 (see text). 



By using two measured frequencies, uji and u n , rather than a single frequency, we obtain a 
measurement of c' independent of the surface tension. Fig. [3] (left panel) shows the values 
of c' that we obtain from the above 'two-frequency' method using the lowest order mode 
n = 2 and a second mode I > 2. We use the n = 2 mode since the corresponding peak 
in the power spectrum is always clearly resolved. We plot the mean of the values obtained 
for each mode / > 2. Error bars (standard error) indicate the variation in the measured d Q 
obtained from different modes I. Filled circles and open squares show data obtained from a 
drop in a nitrogen atmosphere at B Q = 16.5 T, and in air at B = 16.2 T, respectively. The 
broken line on Fig. [3] shows the value of c (r) computed from the potential U shown in Fig. [I] 
(i.e. from the geometry of the solenoid and the current Bq). This line is in reasonably good 
agreement with the experimentally measured values of c , although the data points fall at 
slightly higher values. As an additional check of our method, we plot the dependence of the 
measured c' on I (right hand panel). If the trap-induced frequency enhancement that we 



measure experimentally is accurately given by Eqn. 11 , then the c' values that we determine 



by this method should exhibit no dependence on I. Although there is some scatter in the 
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FIG. 4: Values of the j = (spherically-symmetric) and j = 2, 3 (quadrupole and octopole) coeffi- 
cients c'j in a multipole expansion of the potential well gradient, Eqn.jZJ at Bq = 16.5 T (in nitrogen 
atmosphere). Note, = (see text). The values were determined from the magnetogravitational 
potential (Eqn. [2J shown in Fig. [Tj The Cq line is the same as that shown by the dotted line on 
Fig. i 

data due to experimental error, there is no clear dependence on I. 

V. NON-SPHERICAL POTENTIAL 

In the previous section, we assumed that the potential well was spherically symmetric 
around Lq. However, whilst we have chosen the field Bq to minimize the quadrupole com- 
ponent of the trap, the octopole component remains significant, as can be seen clearly in 
Fig. [TJ Thus it initially appears surprising that we can treat the well as if it were spherical, 
in order to calculate its effect on the vibrations of the droplet's shape. (Had we determined 
an unphysical dependence of our measured c' on / in the previous section, it would also have 
followed that the spherical- well approximation was inadequate.) We now consider the effect 
on the eigenfrequencies of the drop of additional harmonic components c'a of the potential 
trap, and consider why, if the trap has a significant octopole component, the spherical well 
approximation is so effective at reproducing the measured frequencies of the droplet. 

We first discuss how close to spherical it is possible to make the potential trap; that 
is, how small the coefficients c'j, j > can be made. Lorin and Mailfert have considered 
the relationship between the coil geometry (the distribution of the current density) and the 
shape of the potential trap [H]. Although we cannot alter the coil geometry, we can alter 
the coefficients c'j by adjusting the current flowing in the solenoid. The magnitude of the 
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FIG. 5: The three spheres show the computed magnitude of the vector T = — VC/ on the surface of 
three droplets, radius a = 5.0 mm, 7.5 mm and 10.0 mm at Bq = 16.5 T (in nitrogen atmosphere); 
the length and direction of the radiating lines show the direction and magnitude of the (inward 
pointing) T at the surface; T was determined from the magnetogravitational potential (Eqn. [2]) 
shown in Fig. [TJ The variation in |T| over the surface of the droplet is due to the octopole component 
of the potential trap. 

quadrupole component of the radial force \d 2 \ is smallest (for our magnet) at Bq = 16.2 T in 
air, and at B Q = 16.5 T in nitrogen gas; at this field, |c' 2 | <C |c' | and only weakly dependent 
on r. Increasing B beyond this value gives a positive c' 2 , increasing with r. Decreasing 
B below this value gives a negative c' 2 , becoming more negative with increasing r. The 
gradient of the octopole component c' 3 , however, is of the same order as c' and increases as 
r 2 . This component is a feature of the magnetogravitational potential trap generated by a 
solenoid. It can be reduced slightly, by increasing Bq, but only at the expense of significantly 
increasing |c' 2 |. Note that the dipole component of the trap c\ is necessarily zero (hence c[ = 
too), since the net vertical force on the droplet F z must be zero for levitation. Since F z is 
proportional to the difference between the mean surface potentials of its upper and lower 
hemispheres, i.e. F z oc \w\U(a, w)dw — J° 1 \w\U(a,w)dw oc J2j c j /_i wPj(w)dw oc ci, it 
follows that c\ = 0, where w = cos8(= Pi). The harmonic components j > 3, are small 
compared to both c' and c' 3 , in our system. Fig. [4] shows the values of the j = and 
j = 2, 3 (quadrupole and octopole) coefficients c'j at B Q = 16.5 T (for nitrogen atmosphere) 
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computed from the solenoid geometry. The spheres in Fig. [5] show the computed magnitude 
of r = — VC/ at the surface of three droplets, radius a = 5.0 mm, 7.5 mm and 10.0 mm. 
Radiating lines indicate the direction and magnitude of T, which is directed toward the 
interior of the droplet. The variation in |T| over the surface of the droplet is due to the 
octopole component of the potential trap. 

We now consider the effect on the droplet 's eigenfrequencies of adding a harmonic com- 
ponent c'j, j > 2 to the potential trap. The analysis proceeds as above, however, we must 
now include higher order harmonics in the eigenfunction of the shape oscillation 

r = R(6, t) = a + e sin ut 6jPj (cos 9) , (14) 

since / is not, in general, a good eigennumber in a non-spherical potential. In principle, 
we should decompose the shape into spherical harmonics Y™, since the degeneracy in m 
is also lifted in a non-spherical potential (see Fig. [6]). However, our method of inducing 
shape oscillations in the droplet tends to excite only the axisymmetric shapes (i.e. with 
m = 0), since the air jet is aligned along the solenoid axis. For this reason, we derive here 



the frequencies of the m = oscillations only (Eqn. 14), which are sufficient to interpret 
the experimental results. We summarize the treatment of the general case \m\ < I in the 
Appendix. 

The velocity potential is 

(f>(r,6) = -too cos ut ^2bir l l- l a- l+1 Pi{cos9). (15) 

i>i 



The magnetogravitational potential U at the surface of the drop is (see Eqn. 9a) 



U(R) = U(a) - er r (a, 9) sm{ut) ^ hP^cos 9), (16) 



i>i 



where —T r (a, 9) = c' (a) + c'j (a) Pj (cos 9) in this case (see Eqn.[4j). Inserting Eqn. 16, Eqn. 
and Eqn. [15] into Eqn. [61 and equating the time- varying terms, we obtain 
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auj l^ — = 2^ 

i>i i>i l 



(c' (a) + c' j (a)P j ) + ^-(l-l)(l + 2) 



pa 2 



kPi. (17) 



The product P\Pj appearing on the RHS of this equation can be expanded as a sum of 
Legendre polynomials [15J, which, for our purposes, is most conveniently written P\Pj = 

ESj-J| W-?' V + 1 - P]/ 2 ) P P> in which 

Wl '" 7 ' J " A(j + l-s) \2j + 2l-2s + l) [ } 
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for integer s and we define Q(l, j, s) = for half-integer s; A(n) = Ix3x5x...x (2n — l)/n!. 
Equating the coefficients of Pi, we obtain 

co\ = c ^b l + — 3 i(i-m+2)b l 

a pa 6 

+ Cj ^- b P Q(p,j,[j+p-m- (19) 

a V=\3-l\ 

This equation has the form of an eigenvalue problem u 2 bi = Hi\b\ = (H^ + Vi\)b\ (using 
the summation convention, I > 1), where H^°\ representing the first two terms on the RHS 



of Eqn. 19 is a diagonal matrix and V, representing the third term on the RHS is not, in 
general, diagonal. Treating V as a perturbation, the first order correction (cr|)W = 6j Vi\b^ 
due to a harmonic component j can be computed analytically to obtain the eigenfrequencies 
to 2 = a% + ctq + ipff^ + ■ ■ ■■ For a quadrupole harmonic j = 2 we obtain 

/ 2Ui) 4(a) Z 2 (/ + l) 

( ^ } = — (2/ + 3)(2/-l)- (20) 

For an octopole harmonic j = 3, and all odd j, (o'])^ = 0. This explains why our spherical- 
well approximation works so well in predicting the eigenfrequency spectrum of the droplet: 
we have minimized c' 2 by careful adjustment of Bq and we expect, from the above analysis, 
that the effect of the octopole harmonic c' 3 (which is comparable to the spherically symmetric 
component c' ) on the measured eigenfrequencies in our experiment, to be minimal. Fig. [6] 
shows the eigenfrequencies determined by using a computer to solve the eigenvalue problem 
numerically (i.e. beyond first order). The eigenfrequencies of an a = 7.5 mm droplet, for 



various c' 2 and c' 3 , are shown (thick lines), along with the first-order result Eqn. 20 for 
comparison (thin lines). In addition, we plot the eigenfrequencies of modes with m ^ 
(broken lines); the calculation of the frequency of these modes is outlined in the Appendix. 
The lower plot of Fig. [6] shows that the effect of the octopole harmonic c' 3 on the measured 
eigenfrequencies is minimal for all m for \c' 3 \ < 1. For a quadrupole harmonic c' 2 , the 
eigenfrequencies depend strongly on m, as shown in the upper plot of Fig. |6j However, 
since we have reduced \c' 2 \ to smaller than w 0.05 as described above, and since our method 
of exciting the oscillations tends to excite only the m = modes, this does not alter our 
explanation of why the spherical-well approximation works so effectively. 
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VI. DISCUSSION AND CONCLUSION 



Although the forces at the droplet's surface are dominated by the surface tension, we 
can clearly observe the influence of the magnetogravitational potential trap on the eigenfre- 
quencies. We have demonstrated that the effect of the trap on the eigenfrequencies closely 
approximates that of a spherically-symmetric potential well. Although the potential well 
has a significant octopole component in additional to the spherically-symmetric component, 
we have calculated that the effect of the octopole component on the eigenfrequencies is small 
and verified this experimentally. The small difference between the values of c' calculated 
from the solenoid geometry and the measurement of c' from the measured droplet oscil- 
lations, may be due to a small error in the calculated value: the calculation of the field 
profile B(x, z), performed by numerical integration of the Biot-Savart equation, is based on 
a thin-shell approximation of the current density in the solenoid coils. We note that the 
discrepancy cannot be explained by the octopole component of the trap, since this only 
marginally reduces the eigenfrequencies (Fig. [6]), which would only reduce the measured 
value of c' slightly. Experiments have shown that nonlinear effects become significant for 
oscillation amplitudes greater than approximately 0.1a [16]. In our experiments, however, 
the oscillation amplitude is smaller than this, and we do not observe these non-linear effects. 

We have shown how to minimize the quadrupole harmonic of the potential trap, by 



adjusting the solenoid current. This allowed us to use the relatively simple formula Eqn. 11 to 
calculate the effect of the trap on the eigenfrequencies. This opens up the possibility of using 
diamagnetic levitation to accurately determine the surface tension of diamagnetic liquids (for 
example, water, and many water-based and organic solutions): subtracting the contribution 
of the trap to the eigenfrequencies reveals the Rayleigh-frequency spectrum, from which 
the surface tension can be obtained directly. A contactless measurement technique has 
many advantages, for example, the ability to measure highly reactive liquids and to achieve 
significant supercooling of the liquid, as demonstrated in experiments on electromagnetically 
levitated liquid metals [17J. Acoustic levitation and the oscillating drop method can be used 
to measure surface properties of small drops (radius < 2 mm) of water and organic liquids 
(e.g. [IHIIIH]). Suspension of small drops in air-flow is also possible (20]. However, since the 
equilibrium droplet shape is distorted significantly from spherical in both of these techniques, 
it is necessary to make accurate measurements of the equilibrium shape in order to correct 
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for the distortion, which introduces additional significant experimental uncertainty [T8H2T] . 
Drops can be bounced on a solid surface or a vibrating liquid bath, the drop and surface 
being separated by a thin layer of air [22J In this case, the eigenfrequency spectrum 
deviates from the Rayleigh spectrum due to the periodic forcing by the oscillating surface 
[2U [25]. An expression, similar to Rayleigh's, for the eigenfrequencies of vibrating drops 
in continuous point-contact with a surface, has been determined |25j but there remains 
some uncertainty over the spectrum of bouncing drops [25] . Although free-fall can be used 
to obtain the Rayleigh frequencies directly (for example, in a drop-tower [2E]), we have 
shown how diamagnetic levitation could offer an alternative, accurate method of measuring 
these frequencies; the levitated droplets are near-spherical at rest and the eigenfrequency 
spectrum is very close to Rayleigh's, as shown in Fig. [2} The relatively small shift to 
higher frequency due to the potential trap can be obtained using a simple spherical-well 



approximation, Eqn. 11 Diamagnetic levitation can levitate cm-size drops (up to ~ 3 cm- 
diameter in our magnet) [21 IU |27] enabling the droplet volume to be measured easily and 
precisely. Currently, there is interest in the temperature dependence of the surface tension 
of supercooled water [28]. The lowest-temperature achieved so far, which used a contact 
technique, is 245 K [29J; nucleation sites on the container walls trigger freezing before colder 
temperatures can be reached. We propose that diamagnetic levitation and the oscillating 
drop technique could be used to obtain measurements of surface tension at temperatures 
less than 245 K. 

It is interesting to consider an analogue between the oscillations of the levitating droplet 
and that of an object in a gravitational field. We can consider the magnetogravitational 
force on a unit mass of water, T = — V£7, as being an effective gravitational field acting on 
the liquid. Reid has obtained an expression for the eigenfrequencies of a gravitating body 
composed of a solid spherical core, radius R\ and density p\ covered by an inviscid liquid 
'mantle' of radius R 2 , density p 2 [3D] (see also Ref. [3"T]): 



^3"°*<' + Vt(iV-) ^-anyj' (21) 

where r\ = R1/R2, p = i] 3 pi + (1 — rf)p 2 is the mean density and G is the gravitational 



constant. Interestingly, Reid's expression becomes equivalent to Eqn. [TT] if we let p 2 — > 
and .Ri —> 0, since this allows the surface effective gravity T(a) to be equated with a fictitious 
point mass m = (A/3)na 3 p = T(a)a 2 /G ~ 1 x 10 8 kg (approximately the mass of 1 mm 3 of 
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nuclear-density material). In this sense, we can view our magnetically- levitated droplet as 
a 'toy model' of an incompressible gravitating body (possessing surface tension), composed 
of a small-diameter, massive core (the fictitious mass), and a low density outer region (the 
water). Note that this result differs from that of a uniform-density self-gravitating body, 
oq = [(8/3)irGpl(l — 1) / (2Z + 1)] 1//2 [32J . Bastrukov has recently obtained a similar spectrum 
for a self-gravitating liquid droplet with radially- varying density distribution p(r < a) oc 1/r, 
having a singular density at the center [33] . 

By adjusting the current in the magnet (i.e. B ) we can change the quadrupole (c' 2 ) 
component of the trap. Although we have sought to minimize this component in this paper, it 
would be interesting to investigate the effect of this quadrupole component on the oscillation 
frequencies experimentally and compare it with the result we derive above. Unlike the odd 
harmonics, the even order harmonics have a significant effect on the eigenfrequencies to 
first order. Using the analogy with a model 'star', adding a quadrupole component to the 
effective gravity T is equivalent to deforming the shape of the core mass from spherical to 
oblate, perhaps due to rotation, for example. 

This project is supported by a Basic Technology Grant from EPSRC, UK; Grant Nos. 
GR/S83005/01 and EP/G037647/1. 



VII. APPENDIX 



We consider the effect of a trap component c'- on an arbitrary shape oscillation 



r = R(0,$,t) = a + esmutJ2Yl h T Y i m { d ^)- 

i>l m=-l 

The corresponding velocity potential is 

i 



(22) 



(23) 



1>1 m=-l 



Inserting Eqn. [22] and Eqn. [23] into Eqn. [6] and equating the time- varying terms, we obtain 

(24) 



au) 



h m V 



1>1 m=-l 



1>1 m =-l 



T 



U\Aj 



The product Y™Pj = N^e^F^Ff can be expanded as jV/V" 1 * E^-i| Q(°,J; m , 1;p) p ™, 
where the Q in this case are Gaunt coefficients, as defined in Ref. [31], P, m are associated 
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Legendre functions and N™ are the corresponding normalization factors. Equating the 
coefficients of Y™, we obtain 21 + 1 eigenvalue problems, i.e. one for each m, similar to 
Eqn.[T9 
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FIG. 6: Top: Effect on the eigenfrequencies of an a = 7.5 mm water droplet of adding a quadrupole 
harmonic c' 2 to the shape of the magnetogravitational potential well (c' Q = 0.36 ms -2 , all other 
c'a = 0). The thick black lines show the eigenfrequencies (squared) of the m = modes (i.e. the 



eigenvalues u 2 of Eqn. 19) corresponding to I = 1 — 5; the eigenvalues of Eqn. 19 were computed 
using MATLAB. (Note that we have plotted w 2 — o\ for clarity). The thin lines show the first-order 



approximation, Eqn. 20, departing from the numerically-computed frequencies at \d 2 \ > 1. Long- 
dash, medium-dash, short-dash, dot-dash and dotted lines show the frequencies of the \m\ = 1 — 5 
modes, respectively (the treatment of the m/0 modes is summarized in the Appendix). Bottom: 
Effect on the eigenfrequencies of an a = 7.5 mm water droplet of adding a harmonic c' 3 to the shape 
of the well (c' = 0.36 ms -2 , all other c'- = 0). 
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